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Phase- Only  Nulling  as  a  Least-Squares 
Approximation  to  Complex  Weight  Nulling 

1.  INTRODUCTION 

There  is  currently  considerable  interest,  arising  from  the  widespread  use  of 
phased  arrays,  in  phase-only  control  of  array  weights  for  adaptive  cancellation  of 
interference*  ^  and  jr  non-adaptive  shaping  of  antenna  patterns.  J  8  To  develop 
efficient  procedures  of  adaptive  phase-only  cancellation,  it  is  useful  to  have  a  clear 
understanding  of  the  end  products  of  such  adaptive  procedures,  namely,  the  kinds 
of  patterns  with  imposed  nulls  that  can  be  obtained  with  phase-only  weight  perturba¬ 
tions.  Also  of  interest,  because  of  their  potential  utility  in  formulating  adaptive 
procedures,  are  deterministic  methods  of  phase-only  null  synthesis.  It  is  for  these 
reasons  that  we  have  undertaken  a  series  of  related  studies  to  investigate  different 
methods  of  phase-only  null  synthesis. 

The  principal  difference  between  null  synthesis  with  phase -only  control  and 

with  complex  (that  is,  combined  phase  and  amplitude)  control,  is  that  the  equations 

for  imposing  pattern  nulls  with  phase-only  control  are  non-linear,  whereas  the  null 

equations  for  complex  weight  perturbations  are  linear.  The  non-linear  equations 

cannot  be  solved  exactly  analytically  and  hence  alternative  methods  must  be  found 

u 

to  obtain  solutions.  In  our  first  study  we  assumed  that  the  required  phase 
(Received  for  publication  6  May  1982) 

(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 

See  References,  page  21.) 


perturbations  of  the  element  weights  are  small  and  under  this  assumption  linearized 
the  phase -only  null  equations.  The  solution  to  the  linearized  equations  is  a  set  of 
phase  perturbations  that  are  reasonable  approximations  to  the  exact  phase  perturba¬ 
tions  in  certain  situations,  namely  low  sidelobes  and  a  small  number  of  imposed 
nulls.  The  cancellation  effectiveness,  however,  is  not  comparable  to  that  achiev¬ 
able  with  complex  perturbations.  In  our  second  study we  developed  an  iterative 
procedure  for  solving  the  phase- only  null  equations  based  on  repeated  linearization 
of  the  equations.  First,  the  simple  linearization  method  is  used  to  form  a  pattern 
with  lower  power  at  prescribed  locations  than  the  original  pattern.  Then,  the  new 
pattern  is  in  turn  taken  as  the  starting  point  and  the  linearization  procedure  applied 
agai  to  obtain  a  pattern  with  still  lower  values  of  the  power  at  the  prescribed 
locations.  The  procedure  is  repeated  until  nulls  as  deep  as  desired  are  obtained. 
This  iterative  method  is  remarkably  successful  in  low  sidelobe  applications  and  even 
in  applications  to  patterns  with  high  sidelobes  and  a  small  number  of  imposed  null 
locations,  although  it  fails  to  work  as  a  general  method.  When  the  iterative  proce¬ 
dure  does  converge,  the  cancellation  effectiveness  obtained  with  the  phase-only 
perturbations  is  almost  fully  comparable  to  that  obtained  with  complex  perturba¬ 
tions.  The  principal  drawback  of  phase-only  null  synthesis  per  se  (that  is,  with 
near-exact  solutions  to  the  non-linear  null  equations)  is  that  any  nulling  of  the 
original  pattern  is  accompanied  by  a  corresponding  raising  of  the  pattern  at  sym¬ 
metric  locations  with  respect  to  the  mainbeam. 

In  this  report  we  investigate  a  third  method  of  phase-only  null  synthesis:  that 
of  fitting  phase-only  perturbaticas  in  a  least-squares  sense  to  complex  perturbations 
calculated  to  impose  pattern  nulls  it  a  prescribed  set  of  sidelobe  directions.  We 
show  that  such  a  least-squares  fit  is  accomplished  very  simply  by  taking  only  the 
phase  part  o'  the  combined  phase  and  amplitude  perturbations  ar.d  ignoring  the 
amplitude  change  in  the  element  weights.  The  method,  although  extremely  simple 
to  implement,  does  not  in  general  give  cancellation  effectiveness  comparable  to  even 
the  simple  linearized  equation  method,  it  does  however,  give  modest  cancellation 
in  high  sidelobe  and/or  multiple  imposed  null  applications  when  on^  or  both  of  the 
other  two  methods  fail  to  work  at  all. 


2.  ANALYSIS 

We  begin  by  proving  the  central  result  of  this  report:  Let  wQn,  n  =  1,2,  ...  N, 
be  a  set  of  linear  array  weights,  and  let  w  ,  n  =  1,2,  ....  N,  be  a  set  of  weights 
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derived  from  the  original  weights  by  multiplication  by  the  complex  factors  pn 
that  is. 


i*. 


w  =  w  p  e 
n  on 


n  =  1,2,  ... ,  N 

i*r 


Then  the  set  of  weights  w  =  w  e 

pn  on 

tions  that  differs  least  in  a  least-squares  sense  from  the  wr  is 


derived  from  wQn  by  phase-only  perturba- 


w  =  w  e 
pn  on 


n  =  1,2,  ....  N  . 


In  other  words,  the  phase-only  perturbed  weights  that  differ  least  from  the  combined 
phase-and-amplitude  perturbed  weights  are  obtained  by  setting  the  phases  of  the 
phase-only  perturbations  equal  to  the  phases  of  the  combined  phase-and-amplitude 
perturbations. 

The  proof  is  straightforward.  Let  E  be  the  sum  of  the  squares  of  the  differ¬ 
ences  between  the  phase-only  perturbed  weights  and  the  phase-and-amplitude 
perturbed  weights.  Then 


Differentiating  E  with  respect  to  any  particular  one  of  the  6  and  setting  the  deriva¬ 
tive  equal  to  zero 
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whereupon  we  obtain 


(1) 


6  = 
n 


for  any  integer  K.  In  particular  we  can  set  K  =  0  thus  obtaining  the  desired  result. 

This  result  can  be  equivalently  expressed  in  terms  of  patterns  instead  of 
weights,  since  the  sum  of  the  squares  of  the  array  weights  equals  the  mean  square 
of  the  array  pattern,  the  mean  being  taken  over  a  full  period  of  the  variable 
u  -  -y  d  sin  $  (d  -  interelement  spacing,  X  -  wavelength,  9  -  pattern  angle  measured 
from  broadside).  Thus  let  Pq  be  a  linear  array  pattern  and  Pj  a  pattern  derived 
from  Pq  by  perturbations  of  the  amplitudes  and  phases  of  the  array  weights.  Then 
the  pattern  Pg  derived  from  PQ  by  phase-only  perturbations  of  the  array  weights 
that  most  closely  matches  P^  in  a  least  mean  square  sense  is  obtained  by  keeping 
only  the  phases  of  the  weight  perturbations  associated  with  P^  and  neglecting  the 
amplitude  perturbations. 

9 

In  a  previous  report  we  derived  expressions  for  the  minimum  weight 
perturbations,  both  complex  and  small  phase-only,  required  to  impose  nulls  in  a 
given  pattern  at  a  set  of  M  specified  sidelobe  directions.  We  showed  that  the  can¬ 
cellation  pattern  corresponding  to  the  complex  weight  perturbations  could  be  inter¬ 
preted  as  the  sum  of  M  beams  of  the  same  shape,  one  directed  at  each  of  the  null 
locations.  For  small  phase-onlv  perturbations,  the  cancellation  pattern  was 
expressible  as  the  sum  of  M  pairs  of  beams,  one  member  of  each  pair  directed  at 
a  null  location,  and  the  other,  of  opposite  sign,  directed  at  the  symmetric  location 
with  respect  to  the  mainbeam. 

if  we  obtain  a  set  of  phase-only  weight  perturbations  as  a  least-squares  fit  to 
a  given  set  of  complex  perturbations  that  corresponds  to  a  sum  of  M  cancellation 
beams,  one  directed  at  each  of  M  null  locations,  it  is  again  the  case  that  for-  small 
perturbations  the  phase-only  cancellation  pattern  can  be  interpreted  as  the  sum  of 
M  pairs  of  beams  of  opposite  sign.  I'or,  following  Kef.  9,  suppose  that  the  original 
weights  are  given  by 


-id  u 
J  n  s 


n  =  1,2,  ....  N 


(2) 


where  N  is  the  number  of  elements,  the  a  are  the  (real)  amplitudes  of  the  weights, 

^  n_ 

and  we  have  assumed  a  uniform  progressive  phase  shift  =  -y  d  sin  9  s  across  the 
array  with  d  the  interelement  spacing  and  9  the  pattern  angle  measured  from 
broadside.  The  phase  reference  center  is  taken  to  coincide  with  the  center  of  the 
array,  and  the  factor  d^  in  the  exponent  of  fclq.  (2)  is  given  by 

dn  =  1  -  (n  -  1),  n  =  1, 2,  .  . .  ,  N  . 


H 


The  amplitudes  are  assumed  to  be  symmetric  with  respect  to  the  phase  reference 
center.  Let  us  further  suppose  that  the  complex  weight  perturbations,  Awn>  that 
we  are  trying  to  match  in  a  least-squares  sense  by  phase-only  perturbations  are 
given  by 


Aw 


M 

£ 


b  e 
m 


-j  d  u 
n  m 


(3) 


with  the  corresponding  cancellation  pattern 
M  N 


r  h  V  ~  Jdn(u'um) 
m  =  1  ^  n  =  1 


c  e 
n 


In  Eq.  (3),  the  coefficients  cn  that  determine  the  shape  of  the  cancellation  beams 

are  assumed  to  be  real  and  symmetric  with  respect  to  the  phase  reference  center; 

u__  =  ^  d  sin  9  ,  m=l,2,  ....  M  where  0  is  the  pattern  location  of  the  mth  null 

m  A  m  m  r 

and  b  ,  m=  1,  2,  . . . ,  M,  is  a  real  coefficient  giving  the  magnitude  and  sign  of  the 
th 

m  cancellation  beam.  To  obtain  the  phase-only  perturbations,  Aw^,  that  differ 
least  in  a  least-squares  sense  from  the  Awn»  we  first,  following  the  proof  of 
Eq.  (1),  express  the  perturbed  weights,  w  =  w  +  Aw  ,  in  the  form 

i v  n  on  n 

w„  =  w  „  P  e  "  so  that 
n  on  n 


w  =  w  +  Aw 
n  on  n 


(aw  \ 
on  / 

r-  M 

on  E  b  e 

on  L  an  ml  m 

M”»  *is»)] 


-j  d  (u  -  u  ) 
J  n  m  s 


jtan 


'(4:) 


where 


M 


a  =  E  brv,  COS  [  d  ~  u  )1 
n  an  m=  1  m  L  n  m  sj 
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and 

c  M 

a  =  -  -2-  s  *>_  si«  r  <i_  ■  ua>]  * 

n  *n  m  =  1  m  L  n  m  8  J 

The  phase-only  perturbed  weights,  wpn>  differing  least  from  the  combined  phase- 
and-amplitude  perturbed  weights  are  then 


For  small  combined  phase-and-amplitude  perturbations,  both  the  a  ^  and  the 
will  be  small  and  we  can  then  approximate  the  w  by 


w  =  w  <  1  +  jfl  ) 
pn  on  Jrn 


so  that  the  phase-only  weight  perturbations  are 
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The  corresponding  phase-only  cancellation  pattern  is  then 
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J  n 


AFp(U>=  5  AwpnC 
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M  N 

■  b»  £  ‘-i  [ 
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and  so  is,  as  claimed,  the  sum  of  M  pairs  of  beams,  one  member  of  each  pair 

directed  at  a  null  location,  u  =  u  ,  and  the  other  member  of  opposite  sign  directed 

at  the  location  u=u  -2(u  -u),  symmetric  with  respect  to  the  main  beam, 

m  m  s 

This  means  that  if  we  try  to  do  phase-only  nulling  by  a  best  least-squares  fit  of  the 
phase-only  perturbed  weights  to  a  set  of  amplitude-and-phase  perturbed  weights, 
that  give  imposed  nulls  at  the  desired  locations,  any  cancellation  of  the  original 
i-  ftern  will  be  accompanied  by  an  increase  in  the  power  pattern  at  directions 
symmetric  (with  respect  to  the  main  beam)  to  the  null  locations.  No  such  increase 
is  present  in  the  pattern  corresponding  to  the  combined  amplitude-and-phase 
perturbed  weights. 

3.  NUMERICAL  RESULTS  AND  DISCUSSION 

In  this  section  we  present  the  results  of  computations  performed  to  examine  the 

performance  of  phase-only  nulling  as  a  least  squares  approximation  to  complex 

(that  is,  combined  phase-and-amplitude)  nulling.  All  computations  were  performed 

for  an  array  of  41  equispaced  isotropic  elements  with  half  wavelength  interelement 

spacing.  The  results  presented  are  of  two  types.  One  is  the  depth  of  the  null  at 

the  imoosed  null  location(s),  and  the  other  is  the  cancellation  achieved  throughout 

9 

a  sector  of  the  pattern.  Following  our  earlier  report  we  define  the  power  cancella¬ 
tion  ratio,  C,  in  the  sector  A0  =  8^  s  8  s  02>  to  be 

..  sM2 
'  s  [>„<«>] z 

where  F  (0  )  is  the  original  pattern  and  F(0)  is  the  perturbed  pattern. 

In  Table  1  we  show  the  depth  of  null  achieved  with  phase-only  perturbations  ob¬ 
tained  as  a  least-squares  approximation  to  complex  nulling,  and  compare  these  null 
depths  to  those  achieved  with  phase-only  perturbations  obtained  by  solving  the 
linearized  null  equations  (see  Kef.  9).  For  reference  we  also  show  the  null  depths 
obtained  with  an  iterative  method  described  in  Ref.  10  which  attempts  to  find  exact 
solutions  of  the  phase-only  nulling  equations  bv  a  process  of  repeated  linearization. 
The  original  weights  are  for  a  40  dB  Chebvshev  pattern.  Following  Ref.  9,  two  types 
of  complex  nulling  are  considered,  one  in  which  the  sum  of  the  squares  of  the  weight 
perturbations  is  minimized,  and  the  other  in  which  the  sum  of  the  squares  of  the 
weight  perturbations  relative  to  the  original  amplitudes  of  the  weights  is  minimized. 
The  results  presented  are  for  a  series  of  equispaced  | in  sin  (0  )]  imposed  null  loca¬ 
tions.  First,  one  imposed  null  at  15.23°  is  considered.  Then  two  simultaneous 


! 


1  I 


imposed  nulls  at  15.  23°  and  15.  78°  are  considered,  and  so  on  up  to  six  simul¬ 
taneous  imposed  null  locations.  It  is  seen  that  in  general  for  up  to  five  simultaneous 
imposed  nulls,  the  depth  of  null  at  the  imposed  null  locations  obtained  with  phase- 
only  nulling  as  a  least-squares  approximation  is  considerably  shallower  than  that 
obtained  bv  solving  tne  linearized  phase-only  nulling  equations.  Only  for  the  six 
null  case  (when  the  linearized  phase-only  solution  is  almost  completely  ineffective 
because  of  the  large  perturbations  that  result  from  interbeam  interference)  are 
deeper  nulls  obtained  with  the  least-squares  approximated  phase  perturbations.  The 
iterative  linearization  procedure  gives  consistently  far  deeper  nulls  than  either  of 
the  other  two  methods. 

In  Table  2  we  show  the  cancellation  ratios  for  the  sectors  corresponding  to  the 
sets  of  imposed  null  locations  used  in  Table  1.  The  entries  in  the  i**1  row  of  the 
table  are  the  cancellation  ratios  obtained  by  the  methods  listed  across  the  top  of  the 
table  for  the  sector  15.  23  ‘  0  '  (6.-,)  with  nulls  imposed  at  15.  23°  and  at  the  loca¬ 

tions  6  =  (0.,F,  j  -  i.  We  note  that,  for  all  except  the  six  null  case,  the  sector 
cancellation  ratios  are  also  shallower  for  the  least-squares  approximated  phase 
perturbations  than  for  the  perturbations  obtained  as  solutions  to  the  linearized  null 
equations.  This  is  to  be  expected  since  effective  sector  cancellation  depends  in 
large  part  on  deep  nulls  being  formed  at  the  imposed  locations  to  "anchor  the  pattern 
down"  as  it  were.  Since  the  combined  phase-and-amplitude  perturbations  from 
which  the  least-squares  approximated  phase-only  perturbations  are  derived  give 
patterns  with  virtually  perfect  nulls  and  very  effective  sector  cancellation  (see 
Figures  1  and  2)  we  are  forced  to  conclude  that  much  of  the  effectiveness  of  com¬ 
plex  nulling  is  lost  when  we  keep  only  the  phase  part  and  discard  the  amplitude  com¬ 
ponent  of  the  perturbations.  This  loss  of  effectiveness  is  not  an  intrinsic  property 
of  phase-only  nulling  per  se,  as  is  clear  from  the  null  depths  and  sector  cancella¬ 
tion  obtained  with  the  iterative  linearization  procedure.  Rather,  it  is  a  deficiency 
of  the  particular  approximation  method. 

Tables  3  and  4,  similar  to  Tables  1  and  2  respectively,  show  null  depths  and 
cancellation  ratios  for  a  series  of  imposed  null  locations  in  a  20  dB  Chebyshev 
pattern.  The  overall  behavior  of  the  results  obtained  is  similar  to  that  for  the 
40  dB  Chebyshev  pattern  case.  For  a  small  number  of  imposed  nulls,  the  least- 
squares  approximated  phase  perturbations  give  shallower  nulls  and  less  effective 
sector  cancellation  than  do  the  perturbations  obtained  as  solutions  to  the  linearized 
phase-only  null  equations.  As  the  number  of  imposed  nulls  increases,  the  linearized 
phase-only  nerturbation  method  breaks  down  because  large  phase  perturbations  are 
required  to  impose  nulls  and  hence  the  small  phase  perturbation  assumption  used 
to  obtain  the  linearized  phase-only  null  equations  becomes  invalid.  (In  general,  the 
20  dB  Chebyshev  taper  pattern  requires  larger  phase  perturbations  to  impose  nulls 
than  does  the  40  dB  pattern  because  the  pattern  must  be  lowered  further  (that  is. 
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perturbed  more)  to  form  deep  nulls.  )  The  iterative  method,  the  first  iteration  of 
which  coincides  with  the  solution  of  the  linearized  null  equations,  also  breaks  down 
and  further  iterations  fail  to  improve  on  the  initial  iteration.  The  least-squares 
approximated  phase  perturbations,  especially  those  based  on  the  complex  nulling 

method  that  minimizes  zCl^w^f  ,  continues  to  give  some  cancellation  for  five 
or  six  imposed  nulls,  but  the  effectiveness  of  nulling  is  not  impressive  and  is,  in 
general,  less  than  the  effectiveness  of  nulling  shown  for  the  40  dB  case. 


4.  CONCLUSIONS 

In  this  report  we  have  investigated  a  method  of  phase-only  nulling  in  linear 
array  patterns.  The  method  consists  of  obtaining  the  set  of  phase  only  weight 
perturbations  that  differs  least  in  a  least-squares  sense  from  a  set  of  complex 
(that  is,  combined  phase  and  amplitude)  perturbations  calculated  analytically  to 
impose  nulls  at  prescribed  locations  in  a  given  pattern.  We  showed  that  the 
desired  phase  perturbations  are  given  by  taking  the  phase  part  only  of  the  complex 
perturbations  and  ignoring  the  change  in  amplitude  of  the  element  weights. 

Calculations  were  performed  to  examine  the  effectiveness  of  this  nulling 
method  and  to  compare  it  with  two  other  phase-only  methods.  One  of  these  is 
based  on  solving  a  set  of  linearized  equations  obtained  from  the  non-linear  phase- 
only  null  equations  under  the  assumption  of  small  phase  perturbations.  The  other 
attempts  to  find  exact  solutions  of  the  non-linear  null  equations  by  a  process  of 
iterative  linearization.  It  was  found  that  the  least-squares  approximated  phase- 
only  perturbations  did  not  give  deep  nulls  at  the  desired  pattern  locations  and,  for 
a  small  number  of  imposed  nulls,  was  less  effective  than  the  linearized  phase-only 
perturbation  method,  which  in  turn  is  usually  considerably  less  effective  than  the 
iterative  procedure.  However,  the  least-squares  approximated  phase  method  did 
provide  some  cancellation  in  situations  where  one  or  both  of  the  other  two  methods 
failed  to  work  at  all;  namely,  patterns  with  high  sidelobes  and/or  many  imposed 
null  locations. 
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Table  1.  Depth  of  Null  (dB)  at  Imposed  Null  Directions  0  40  dB  Chebyshev  Amplitude  Taper 


Table  3.  Depth  of  Null  (dB)  at  Imposed  Null  Directions,  0  20  dB  Chebyshev  Amplitude  Taper 
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Figure  la.  Perturbed  Patterns  Obtained  by  Complex  Weight 

Perturbations  ( - )  and  Least -squares  Approximated 

Phase-only  Perturbations  { - )  With  One  Imposed  Null  at 

15.23°.  Range  of  0:  [-90°,  90°] 


Figure  lb.  Perturbed  Patterns  Obtained  by  Complex  Weight 

Perturbations  ( - )  and  Least-squares  Approximated 

Phase-only  Perturbations  ( - )  With  One  Imposed  Null  at 

15.23°.  Range  of  6:  [10°,  30°) 
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Figure  2a.  Perturbed  Patterns  Obtained  by  Complex  Weight 

Perturbations  ( - )  and  Least-squares  Approximated 

Phase -only  Perturbations  ( - )  With  Six  Imposed  Nulls  at 

15.23°,  15.78°,  16.33°,  16.88°,  17.44°,  and  18.00°. 

Range  of  0:  [-90°,  90°) 
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Figure  2b.  Perturbed  Patterns  Obtained  by  Complex  Weight 

Perturbations  ( - )  and  Least-squares  Approximated 

Phase -only  Perturbations  ( - )  With  Six  Imposed  Nulls  at 

15.23°,  15.78°,  16.33°,  16.88°,  17.44°,  and  18.00°. 

Range  of  9:  [10°,  30°) 
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